Parallel Streams of Direct Corticogeniculate Feedback from Mid-level Extrastriate Cortex in the Macaque Monkey.

First-order thalamic nuclei receive feedforward signals from peripheral receptors and relay these signals to primary sensory cortex. Primary sensory cortex, in turn, provides reciprocal feedback to first-order thalamus. Because the vast majority of sensory thalamocortical inputs target primary sensory cortex, their complementary corticothalamic neurons are assumed to be similarly restricted to primary sensory cortex. We upend this assumption by characterizing morphologically diverse neurons in multiple mid-level visual cortical areas of the primate (Macaca mulatta) brain that provide direct feedback to the primary visual thalamus, the dorsal lateral geniculate nucleus (LGN). Although the majority of geniculocortical neurons project to primary visual cortex (V1), a minority, located mainly in the koniocellular LGN layers, provide direct input to extrastriate visual cortex. These "V1-bypassing" projections may be implicated in blindsight. We hypothesized that geniculocortical inputs directly targeting extrastriate cortex should be complemented by reciprocal corticogeniculate circuits. Using virus-mediated circuit tracing, we discovered corticogeniculate neurons throughout three mid-level extrastriate areas: MT, MST, and V4. Quantitative morphological analyses revealed nonuniform distributions of unique cell types across areas. Many extrastriate corticogeniculate neurons had spiny stellate morphology, suggesting possible targeting of koniocellular LGN layers. Importantly though, multiple morphological types were observed across areas. Such morphological diversity could suggest parallel streams of V1-bypassing corticogeniculate feedback at multiple stages of the visual processing hierarchy. Furthermore, the presence of corticogeniculate neurons across visual cortex necessitates a reevaluation of the LGN as a hub for visual information rather than a simple relay.

High-Resolution Laminar Identification in Macaque Primary Visual Cortex Using Neuropixels Probes.

Laminar electrode arrays allow simultaneous recording of activity of many cortical neurons and assignment to correct layers using current source density (CSD) analyses. Electrode arrays with 100-micron contact spacing can estimate borders between layer 4 versus superficial or deep layers, but in macaque primary visual cortex (V1) there are far more layers, such as 4A which is only 50-100 microns thick. Neuropixels electrode arrays have 20-micron spacing, and thus could potentially discern thinner layers and more precisely identify laminar borders. Here we show that CSD signals lack the spatial resolution required to take advantage of high density Neuropixels arrays and describe the development of approaches based on higher resolution electrical signals and analyses, including spike waveforms and spatial spread, unit density, high-frequency action potential (AP) power spectrum, temporal power change, and coherence spectrum, that afford far higher resolution of laminar distinctions, including the ability to precisely detect the borders of even the thinnest layers of V1.

Parallel pathways carrying direction-and orientation-selective retinal signals to layer 4 of the mouse visual cortex.

Parallel visual pathways from the retina to the primary visual cortex (V1) via the lateral geniculate nucleus are common to many mammalian species, including mice, carnivores, and primates. However, it remains unclear which visual features present in both retina and V1 may be inherited from parallel pathways versus extracted by V1 circuits in the mouse. Here, using calcium imaging and rabies circuit tracing, we explore the relationships between tuning of layer 4 (L4) V1 neurons and their retinal ganglion cell (RGC) inputs. We find that subpopulations of L4 V1 neurons differ in their tuning for direction, orientation, spatial frequency, temporal frequency, and speed. Furthermore, we find that direction-tuned L4 V1 neurons receive input from direction-selective RGCs, whereas orientation-tuned L4 V1 neurons receive input from orientation-selective RGCs. These results suggest that direction and orientation tuning of V1 neurons may be partly inherited from parallel pathways originating in the retina.

Postsynaptic cell type and synaptic distance do not determine efficiency of monosynaptic rabies virus spread measured at synaptic resolution.

Retrograde monosynaptic tracing using glycoprotein-deleted rabies virus is an important component of the toolkit for investigation of neural circuit structure and connectivity. It allows for the identification of first-order presynaptic connections to cell populations of interest across both the central and peripheral nervous system, helping to decipher the complex connectivity patterns of neural networks that give rise to brain function. Despite its utility, the factors that influence the probability of transsynaptic rabies spread are not well understood. While it is well established that expression levels of rabies glycoprotein used to trans-complement G-deleted rabies can result in large changes in numbers of inputs labeled per starter cell (convergence index [CI]), it is not known how typical values of CI relate to the proportions of synaptic contacts or input neurons labeled. And it is not known whether inputs to different cell types, or synaptic contacts that are more proximal or distal to the cell body, are labeled with different probabilities. Here, we use a new rabies virus construct that allows for the simultaneous labeling of pre- and postsynaptic specializations to quantify the proportion of synaptic contacts labeled in mouse primary visual cortex. We demonstrate that with typical conditions about 40% of first-order presynaptic excitatory synapses to cortical excitatory and inhibitory neurons are labeled. We show that using matched tracing conditions there are similar proportions of labeled contacts onto L4 excitatory pyramidal, somatostatin (Sst) inhibitory, and vasoactive intestinal peptide (Vip) starter cell types. Furthermore, we find no difference in the proportions of labeled excitatory contacts onto postsynaptic sites at different subcellular locations.

Conserved and divergent gene regulatory programs of the mammalian neocortex.

Divergence of cis-regulatory elements drives species-specific traits1, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains unclear. Here we investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset and mouse using single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome and chromosomal conformation profiles from a total of over 200,000 cells. From these data, we show evidence that divergence of transcription factor expression corresponds to species-specific epigenome landscapes. We find that conserved and divergent gene regulatory features are reflected in the evolution of the three-dimensional genome. Transposable elements contribute to nearly 80% of the human-specific candidate cis-regulatory elements in cortical cells. Through machine learning, we develop sequence-based predictors of candidate cis-regulatory elements in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Finally, we show that epigenetic conservation combined with sequence similarity helps to uncover functional cis-regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.

Asymmetric cortical projections to striatal direct and indirect pathways distinctly control actions.

The striatal direct and indirect pathways constitute the core for basal ganglia function in action control. Although both striatal D1- and D2-spiny projection neurons (SPNs) receive excitatory inputs from the cerebral cortex, whether or not they share inputs from the same cortical neurons, and how pathway-specific corticostriatal projections control behavior remain largely unknown. Here using a new G-deleted rabies system in mice, we found that more than two-thirds of excitatory inputs to D2-SPNs also target D1-SPNs, while only one-third do so vice versa. Optogenetic stimulation of striatal D1- vs. D2-SPN-projecting cortical neurons differently regulate locomotion, reinforcement learning and sequence behavior, implying the functional dichotomy of pathway-specific corticostriatal subcircuits. These results reveal the partially segregated yet asymmetrically overlapping cortical projections on striatal D1- vs. D2-SPNs, and that the pathway-specific corticostriatal subcircuits distinctly control behavior. It has important implications in a wide range of neurological and psychiatric diseases affecting cortico-basal ganglia circuitry.

Parallel pathways carrying direction and orientation selective retinal signals to layer 4 of mouse visual cortex.

Parallel functional and anatomical visual pathways from the retina to primary visual cortex (V1) via the lateral geniculate nucleus (LGN) are common to many mammalian species, including mice, carnivores and primates. However, the much larger number of retinal ganglion cell (RGC) types that project to the LGN, as well as the more limited lamination of both the LGN and the thalamocortical-recipient layer 4 (L4) in mice, leaves considerable uncertainty about which visual features present in both retina and V1 might be inherited from parallel pathways versus extracted by V1 circuits in the mouse visual system. Here, we explored the relationships between functional properties of L4 V1 neurons and their RGC inputs by taking advantage of two Cre-expressing mouse lines - Nr5a1-Cre and Scnn1a-Tg3-Cre - that each label functionally and anatomically distinct populations of L4 neurons. Visual tuning properties of L4 V1 neurons were evaluated using Cre-dependent expression of GCaMP6s followed by 2-photon calcium imaging. RGCs providing input to these neurons (via LGN) were labeled and characterized using Cre-dependent trans-synaptic retrograde labeling with G-deleted rabies virus. We find significant differences in the tuning of Nr5a1-Cre versus Scnn1a-Tg3-Cre neurons for direction, orientation, spatial frequency, temporal frequency, and speed. Strikingly, a subset of the RGCs had tuning properties that matched the direction and orientation tuning properties of the L4 V1 neurons to which they provided input. Altogether, these results suggest that direction and orientation tuning of V1 neurons may be at least partly inherited from parallel pathways originating in the retina.

Transcriptomic cytoarchitecture reveals principles of human neocortex organization.

Variation in cytoarchitecture is the basis for the histological definition of cortical areas. We used single cell transcriptomics and performed cellular characterization of the human cortex to better understand cortical areal specialization. Single-nucleus RNA-sequencing of 8 areas spanning cortical structural variation showed a highly consistent cellular makeup for 24 cell subclasses. However, proportions of excitatory neuron subclasses varied substantially, likely reflecting differences in connectivity across primary sensorimotor and association cortices. Laminar organization of astrocytes and oligodendrocytes also differed across areas. Primary visual cortex showed characteristic organization with major changes in the excitatory to inhibitory neuron ratio, expansion of layer 4 excitatory neurons, and specialized inhibitory neurons. These results lay the groundwork for a refined cellular and molecular characterization of human cortical cytoarchitecture and areal specialization.

Comparative single cell epigenomic analysis of gene regulatory programs in the rodent and primate neocortex.

Sequence divergence of cis- regulatory elements drives species-specific traits, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains to be elucidated. We investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset, and mouse with single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome, and chromosomal conformation profiles from a total of over 180,000 cells. For each modality, we determined species-specific, divergent, and conserved gene expression and epigenetic features at multiple levels. We find that cell type-specific gene expression evolves more rapidly than broadly expressed genes and that epigenetic status at distal candidate cis -regulatory elements (cCREs) evolves faster than promoters. Strikingly, transposable elements (TEs) contribute to nearly 80% of the human-specific cCREs in cortical cells. Through machine learning, we develop sequence-based predictors of cCREs in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Lastly, we show that epigenetic conservation combined with sequence similarity helps uncover functional cis -regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.

Monosynaptic restriction of the anterograde herpes simplex virus strain H129 for neural circuit tracing.

Identification of synaptic partners is a fundamental task for systems neuroscience. To date, few reliable techniques exist for whole brain labeling of downstream synaptic partners in a cell-type-dependent and monosynaptic manner. Herein, we describe a novel monosynaptic anterograde tracing system based on the deletion of the gene UL6 from the genome of a cre-dependent version of the anterograde Herpes Simplex Virus 1 strain H129. Given that this knockout blocks viral genome packaging and thus viral spread, we reasoned that co-infection of a HSV H129 ΔUL6 virus with a recombinant adeno-associated virus expressing UL6 in a cre-dependent manner would result in monosynaptic spread from target cre-expressing neuronal populations. Application of this system to five nonreciprocal neural circuits resulted in labeling of neurons in expected projection areas. While some caveats may preclude certain applications, this system provides a reliable method to label postsynaptic partners in a brain-wide fashion.

Cone opponent functional domains in primary visual cortex combine signals for color appearance mechanisms.

Studies of color perception have led to mechanistic models of how cone-opponent signals from retinal ganglion cells are integrated to generate color appearance. But it is unknown how this hypothesized integration occurs in the brain. Here we show that cone-opponent signals transmitted from retina to primary visual cortex (V1) are integrated through highly organized circuits within V1 to implement the color opponent interactions required for color appearance. Combining intrinsic signal optical imaging (ISI) and 2-photon calcium imaging (2PCI) at single cell resolution, we demonstrate cone-opponent functional domains (COFDs) that combine L/M cone-opponent and S/L + M cone-opponent signals following the rules predicted from psychophysical studies of color perception. These give rise to an orderly organization of hue preferences of the neurons within the COFDs and the generation of hue "pinwheels". Thus, spatially organized neural circuits mediate an orderly transition from cone-opponency to color appearance that begins in V1.

Diversity in spatial frequency, temporal frequency, and speed tuning across mouse visual cortical areas and layers.

The mouse visual system consists of several visual cortical areas thought to be specialized for different visual features and/or tasks. Previous studies have revealed differences between primary visual cortex (V1) and other higher visual areas, namely, anterolateral (AL) and posteromedial (PM), and their tuning preferences for spatial and temporal frequency. However, these differences have primarily been characterized using methods that are biased toward superficial layers of cortex, such as two-photon calcium imaging. Fewer studies have investigated cell types in deeper layers of these areas and their tuning preferences. Because superficial versus deep-layer neurons and different types of deep-layer neurons are known to have different feedforward and feedback inputs and outputs, comparing the tuning preferences of these groups is important for understanding cortical visual information processing. In this study, we used extracellular electrophysiology and two-photon calcium imaging targeted toward two different layer 5 cell classes to characterize their tuning properties in V1, AL, and PM. We find that deep-layer neurons, similar to superficial layer neurons, are also specialized for different spatial and temporal frequencies, with the strongest differences between AL and V1, and AL and PM, but not V1 and PM. However, we note that the deep-layer neuron populations preferred a larger range of SFs and TFs compared to previous studies. We also find that extratelencephalically projecting layer 5 neurons are more direction selective than intratelencephalically projecting layer 5 neurons.

Targeting thalamic circuits rescues motor and mood deficits in PD mice.

Although bradykinesia, tremor and rigidity are the hallmark motor defects in patients with Parkinson's disease (PD), patients also experience motor learning impairments and non-motor symptoms such as depression1. The neural circuit basis for these different symptoms of PD are not well understood. Although current treatments are effective for locomotion deficits in PD2,3, therapeutic strategies targeting motor learning deficits and non-motor symptoms are lacking4-6. Here we found that distinct parafascicular (PF) thalamic subpopulations project to caudate putamen (CPu), subthalamic nucleus (STN) and nucleus accumbens (NAc). Whereas PF→CPu and PF→STN circuits are critical for locomotion and motor learning, respectively, inhibition of the PF→NAc circuit induced a depression-like state. Whereas chemogenetically manipulating CPu-projecting PF neurons led to a long-term restoration of locomotion, optogenetic long-term potentiation (LTP) at PF→STN synapses restored motor learning behaviour in an acute mouse model of PD. Furthermore, activation of NAc-projecting PF neurons rescued depression-like phenotypes. Further, we identified nicotinic acetylcholine receptors capable of modulating PF circuits to rescue different PD phenotypes. Thus, targeting PF thalamic circuits may be an effective strategy for treating motor and non-motor deficits in PD.

Single-cell transcriptomic classification of rabies-infected cortical neurons.

Cortical circuit tracing using modified rabies virus can identify input neurons making direct monosynaptic connections onto neurons of interest. However, challenges remain in our ability to establish the cell type identity of rabies-labeled input neurons. While transcriptomics may offer an avenue to characterize inputs, the extent of rabies-induced transcriptional changes in distinct neuronal cell types remains unclear, and whether these changes preclude characterization of rabies-infected neurons according to established transcriptomic cell types is unknown. We used single-nucleus RNA sequencing to survey the gene expression profiles of rabies-infected neurons and assessed their correspondence with established transcriptomic cell types. We demonstrated that when using transcriptome-wide RNA profiles, rabies-infected cortical neurons can be transcriptomically characterized despite global and cell-type-specific rabies-induced transcriptional changes. Notably, we found differential modulation of neuronal marker gene expression, suggesting that caution should be taken when attempting to characterize rabies-infected cells with single genes or small gene sets.

Secondary auditory cortex mediates a sensorimotor mechanism for action timing.

The ability to accurately determine when to perform an action is a fundamental brain function and vital to adaptive behavior. The behavioral mechanism and neural circuit for action timing, however, remain largely unknown. Using a new, self-paced action timing task in mice, we found that deprivation of auditory, but not somatosensory or visual input, disrupts learned action timing. The hearing effect was dependent on the auditory feedback derived from the animal's own actions, rather than passive environmental cues. Neuronal activity in the secondary auditory cortex was found to be both correlated with and necessary for the proper execution of learned action timing. Closed-loop, action-dependent optogenetic stimulation of the specific task-related neuronal population within the secondary auditory cortex rescued the key features of learned action timing under auditory deprivation. These results unveil a previously underappreciated sensorimotor mechanism in which the secondary auditory cortex transduces self-generated audiomotor feedback to control action timing.

DNA methylation atlas of the mouse brain at single-cell resolution.

Mammalian brain cells show remarkable diversity in gene expression, anatomy and function, yet the regulatory DNA landscape underlying this extensive heterogeneity is poorly understood. Here we carry out a comprehensive assessment of the epigenomes of mouse brain cell types by applying single-nucleus DNA methylation sequencing1,2 to profile 103,982 nuclei (including 95,815 neurons and 8,167 non-neuronal cells) from 45 regions of the mouse cortex, hippocampus, striatum, pallidum and olfactory areas. We identified 161 cell clusters with distinct spatial locations and projection targets. We constructed taxonomies of these epigenetic types, annotated with signature genes, regulatory elements and transcription factors. These features indicate the potential regulatory landscape supporting the assignment of putative cell types and reveal repetitive usage of regulators in excitatory and inhibitory cells for determining subtypes. The DNA methylation landscape of excitatory neurons in the cortex and hippocampus varied continuously along spatial gradients. Using this deep dataset, we constructed an artificial neural network model that precisely predicts single neuron cell-type identity and brain area spatial location. Integration of high-resolution DNA methylomes with single-nucleus chromatin accessibility data3 enabled prediction of high-confidence enhancer-gene interactions for all identified cell types, which were subsequently validated by cell-type-specific chromatin conformation capture experiments4. By combining multi-omic datasets (DNA methylation, chromatin contacts, and open chromatin) from single nuclei and annotating the regulatory genome of hundreds of cell types in the mouse brain, our DNA methylation atlas establishes the epigenetic basis for neuronal diversity and spatial organization throughout the mouse cerebrum.

Epigenomic diversity of cortical projection neurons in the mouse brain.

Neuronal cell types are classically defined by their molecular properties, anatomy and functions. Although recent advances in single-cell genomics have led to high-resolution molecular characterization of cell type diversity in the brain1, neuronal cell types are often studied out of the context of their anatomical properties. To improve our understanding of the relationship between molecular and anatomical features that define cortical neurons, here we combined retrograde labelling with single-nucleus DNA methylation sequencing to link neural epigenomic properties to projections. We examined 11,827 single neocortical neurons from 63 cortico-cortical and cortico-subcortical long-distance projections. Our results showed unique epigenetic signatures of projection neurons that correspond to their laminar and regional location and projection patterns. On the basis of their epigenomes, intra-telencephalic cells that project to different cortical targets could be further distinguished, and some layer 5 neurons that project to extra-telencephalic targets (L5 ET) formed separate clusters that aligned with their axonal projections. Such separation varied between cortical areas, which suggests that there are area-specific differences in L5 ET subtypes, which were further validated by anatomical studies. Notably, a population of cortico-cortical projection neurons clustered with L5 ET rather than intra-telencephalic neurons, which suggests that a population of L5 ET cortical neurons projects to both targets. We verified the existence of these neurons by dual retrograde labelling and anterograde tracing of cortico-cortical projection neurons, which revealed axon terminals in extra-telencephalic targets including the thalamus, superior colliculus and pons. These findings highlight the power of single-cell epigenomic approaches to connect the molecular properties of neurons with their anatomical and projection properties.

Distinct "driving" versus "modulatory" influences of different visual corticothalamic pathways.

Higher-order (HO) thalamic nuclei interact extensively and reciprocally with the cerebral cortex. These corticothalamic (CT) interactions are thought to be important for sensation and perception, attention, and many other important brain functions. CT projections to HO thalamic nuclei, such as the visual pulvinar, originate from two different excitatory populations in cortical layers 5 and 6, whereas first-order nuclei (such as the dorsolateral geniculate nucleus; dLGN) only receive layer 6 CT input. It has been proposed that these layer 5 and layer 6 CT pathways have different functional influences on the HO thalamus, but this has never been directly tested. By optogenetically inactivating different CT populations in the primary visual cortex (V1) and recording single-unit activity from V1, dLGN, and pulvinar of awake mice, we demonstrate that layer 5, but not layer 6, CT projections drive visual responses in the pulvinar, even while both pathways provide retinotopic, baseline excitation to their thalamic targets. Inactivating the superior colliculus also suppressed visual responses in the same subregion of the pulvinar, demonstrating that cortical layer 5 and subcortical inputs both contribute to HO visual thalamic activity-even at the level of putative single neurons. Altogether, these results indicate a functional division of "driver" and "modulator" CT pathways from V1 to the visual thalamus in vivo.

Application of Recombinant Rabies Virus to Xenopus Tadpole Brain.

The Xenopus laevis experimental system has provided significant insight into the development and plasticity of neural circuits. Xenopus neuroscience research would be enhanced by additional tools to study neural circuit structure and function. Rabies viruses are powerful tools to label and manipulate neural circuits and have been widely used to study mesoscale connectomics. Whether rabies virus can be used to transduce neurons and express transgenes in Xenopus has not been systematically investigated. Glycoprotein-deleted rabies virus transduces neurons at the axon terminal and retrogradely labels their cell bodies. We show that glycoprotein-deleted rabies virus infects local and projection neurons in the Xenopus tadpole when directly injected into brain tissue. Pseudotyping glycoprotein-deleted rabies with EnvA restricts infection to cells with exogenous expression of the EnvA receptor, TVA. EnvA pseudotyped virus specifically infects tadpole neurons with promoter-driven expression of TVA, demonstrating its utility to label targeted neuronal populations. Neuronal cell types are defined by a combination of features including anatomical location, expression of genetic markers, axon projection sites, morphology, and physiological properties. We show that driving TVA expression in one hemisphere and injecting EnvA pseudotyped virus into the contralateral hemisphere, retrogradely labels neurons defined by cell body location and axon projection site. Using this approach, rabies can be used to identify cell types in Xenopus brain and simultaneously to express transgenes which enable monitoring or manipulation of neuronal activity. This makes rabies a valuable tool to study the structure and function of neural circuits in Xenopus.Significance StatementStudies in Xenopus have contributed a great deal to our understanding of brain circuit development and plasticity, regeneration, and hormonal regulation of behavior and metamorphosis. Here, we show that recombinant rabies virus transduces neurons in the Xenopus tadpole, enlarging the toolbox that can be applied to studying Xenopus brain. Rabies can be used for retrograde labeling and expression of a broad range of transgenes including fluorescent proteins for anatomical tracing and studying neuronal morphology, voltage or calcium indicators to visualize neuronal activity, and photo- or chemosensitive channels to control neuronal activity. The versatility of these tools enables diverse experiments to analyze and manipulate Xenopus brain structure and function, including mesoscale connectivity.